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Abstract 

The NuTeV Collaboration has observed an excess in their di-muon channel, possibly 
corresponding to a long-lived neutral particle with only weak interactions and which 
decays to muon pairs. We show that this can not be explained by pair production 
of neutralinos in the target followed by their decay far downstream in the detector 
via a LLE R-parity violating (RPV) operator, as suggested in the literature. In 
the parameter region allowed by LEP the event rate is far too small. We propose 
instead a new neutralino production method via U-mesons, which can fully explain 
the observation. This is analogous to neutrino production via 7r-mesons. This 
model can be completely tested and thus also possibly excluded with NOMAD 
data. If it is excluded, the NuTeV observation is most likely not due to physics 
beyond the Standard Model. Our model can also be tested at the current and 
future i?-factories. This opens up a new way of testing a long-lived neutralino LSP 
at fixed-target experiments and thus a possibility to close the gap between collider 
and cosmological tests of R-parity violation. We also discuss a possible explanation 
in terms of a neutral heavy lepton mixing with the Standard Model neutrinos. The 
flavour structure of the observation can be accounted for but the production rate is 
far too low. 



1 Introduction 

In supersymmetry with broken R-parity the MSSM superpotential is extended 

by 

Wjtp = XijkCabLiLjEk + X^eabL^QjDk + A-^e^f/- D D k + ^e a6 L"if*. (1) 

Here L, Q (E, U, D) are the lepton and quark doublet (singlet) left-handed chiral su- 
perfields, respectively. A, A', A" are dimensionless coupling constants, i,j, k = 1,2,3 are 
generation indices. a,b = 1,2 and a, j3, 7 = 1, 2, 3 are SU (2)l and SU(3) C gauge indices, 
respectively. The main phenomenological changes to the MSSM are that the lightest 
supersymmetric particle (LSP) is no longer stable and supersymmetric particles can be 



produced singly at colliders. Through resonance production the couplings (A, A', A") can 
be probed down to about 10~ 3 before the production cross section becomes too small 
If we consider MSSM supersymmetric pair production with a neutralino LSP then we can 
typically probe couplings down to 10~ 5 or 1CT 6 [p|-|i"2|j. For smaller couplings the LSP 
decays outside the detector and we retrieve the MSSM signatures at colliders. Cosmo- 
logically one can exclude lifetimes for the LSP between 1 s < r x o < 10 17 yr which 
corresponds to couplings 10~ 22 < (A, A', A") < 10 -10 . This leaves a gap in experimental 
sensitivity to the R-parity violating couplings^ 10 -10 < (A, A', A") < 1CT 6 [0. Fixed-target 
experiments with remote detectors can probe significantly longer lifetimes than collider 
experiments and are thus an ideal environment for closing this gap in sensitivity |[L4|| . 

The NuTeV Collaboration has searched for long-lived neutral particles (N ) with mass 
Mjyo > 2.2 GeV and small interaction rates with ordinary matter p~5HT7| . They look for 
the decay of the neutral particles in a detector which is 1.4 km downstream from the 
production point. They observe 3 /i/z events where they only expect to see a background 
of 0.069 ± 0.010 events. The probability that this is a fluctuation of this specific channel 
is about 8 • 10~ 5 , which corresponds to about 4.6 a. The probability for a fluctuation of 
this magnitude into any of the di-lepton channels is about 3er. 

The NuTeV experiment considered in detail the possibility that this discrepancy 
is due to a N° which decays into a three-body final state. In Ref. |16| several kine- 



matic distributions of the di-muon events were checked against the hypothesis of a iV° 
with mass 5 GeV: the transverse mass-, invariant di-muon mass- and the missing p?- 
distributions all agree well with the N° hypothesis. The distribution in the energy asym- 
metry Ae = | Ei — E 2 \/(Ei + E 2 ) of the three events (Ei and E 2 are the two observed 
muon energies in each event) shows a low probability for the iV° hypothesis. Thus three 
out of four distributions work very well and as NuTeV does, we consider it worthwhile 
to investigate whether this observation could be due to new physics. It is the purpose 
of this letter to consider two possible models which could explain the observation: (i) a 
light neutralino which decays via R-parity violation and (ii) a neutral heavy lepton (NHL) 
mixing with the Standard Model neutrinos. 

A search for R-parity violating neutralino decays at NuTeV has been proposed in 
Ref. jnj and the couplings A122 and A133 were discussed. In Ref. |i~6l , |i~7| NuTeV themselves 
mention the possibility of R-parity violating neutralino decays as a solution to the observed 



discrepancy, without looking at any specific couplings. In Ref. |T5[ NuTeV searched for 
the neutralino of a very specific model ||18|| . This neutralino was very light and decayed 
via L\L 2 E\ or L\Lj,E\ to an ee final state. Certain supersymmetric parameter ranges 
were excluded assuming neutralino pair-production. 

Here we show that the simple scenarios discussed in the literature can not lead to an 
excess at NuTeV, since the decisive supersymmetric parameter range to get a significant 
neutralino production cross section has been excluded by LEP. We propose instead the 
production of light neutralinos via 5-mesons which could give a measurable excess. We 
briefly present the two possible models and then discuss them quantitatively. 

In section |5] we show that the production rate for neutral heavy leptons is also too low 
and does not lead to a viable explanation. 



1 These coupling values have been determined for a photino LSP of M x o = O(50) GeV and scalar 
fermion masses of Mr = O(100 GeV). 




Figure 1: Neutralino decays through the R-parity violating coupling A 2 32- Diagrams (a-c) 
give rise to di-muon events while the diagrams (d-f) to tau-muon ones. The index a— 1,2 
denotes the mass eigenstate of the slepton. 

2 The Rp Violating Model 

The heavy neutral particle we consider is the lightest neutralino x?, which we also assume 
to be the LSP. In the notation of 0, the neutralino decays as Xi ~^ P#p> where 0^ p is 
the dominant R-parity violating operator. Only two operators give a di-muon signature: 
^2i2tabL%L\E'i, % = 1,3. For i = 1 the neutralino will decay with equal probability to 
e/iu and \x\xv. No e/i-events are observed, we therefore propose one dominant R-parity 
violating operator: 

Otfj, = tab^L^LlE^. (2) 
For later reference we quote the experimental bound on this operator [19] 



A 232 < 0.070 ( m ** ) , (2a). (3) 
3 VIOOGeV/ V ; V ; 

The operator in Eq.(§) corresponds to the two neutralino decay modes (Fig.l) 

x o _ | (4) 

as well as their complex conjugate, since the neutralino is a Majorana spinor. We shall 
show below that for a light neutralino the t\i decays are sufficiently phase space suppressed 
to give an expectation below one event. For the light neutralino production we shall 
consider two possibilities: 

1. Pair production of the neutralinos |2(J which proceeds via (a) s-channel Z° boson 
exchange and (b) t-channel squark exchange. 

2. Single neutralino production in the decay of bottom hadrons. The bottom hadrons 
are formed following the production of a bb pair. These hadrons can then decay via 
the R-parity couplings X' m , (i = 1, 2, 3). We will only consider the decays of the 
and B + via R-parity violation (Fig. 2) 
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Figure 2: Neutralino production in 5-meson decays : (a-c) — > ViXi, an d (d-f) -B" 



This mechanism allows one to produce light neutralinos via a strong interaction 
process and is analogous to the production of neutrino beams via 7r's and K's 
(and D's). A related mechanism was discussed in the context of the Karmen time 
anomaly [ R| 21 1 . 



For later reference we present the experimental bounds on the A- 13 at 2 a P,|19 
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3 Quantitative Analysis 

As discussed by the NuTeV experiment, the mass of the N° is roughly 5 GeV. The con- 
straints on a very light neutralino were discussed in detail in Ref. |18j. We expect them 
to mainly carry over to the present mass region |H. In order to get a M x? = 0(5 GeV) 
neutralino and avoid the LEP bounds we must consider the case, where the electroweak 
gaugino masses Mi, M^ are independent parameters. In Fig. ||] we show the MSSM pa- 
rameter space which corresponds to M x o = (5 ± 0.5) GeV for two values of tan/3 and 
sgn jj,. The composition of the neutralino is more than 99% bino, provided the lightest 
chargino mass is greater than 100 GeV. 

The dominant bino-nature of the LSP has immediate implications for pair production 
of neutralinos. The bino does not couple to the Z° boson and thus the s-channel pair- 
production of the bino is negligible. This only leaves the t-channel production which is 
proportional to M^ 4 and thus strongly suppressed. We shall quantify this below. 

In both cases neutralino production is followed by the decay. The matrix elements 
for the decay via Jji p were given in [p], |23|. As the neutralino in our model will be 
much lighter than the sleptons (Mg > 90 GeV from LEP) it is sufficient to neglect the 
momentum flow through the slepton propagators. For a purely bino neutralino in this 
limit the spin averaged matrix element is given by 
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Figure 3: Solutions in (M 1; M 2 , fi, tan/5) giving 4.5 GeV < M^o < 5.5 GeV in the cross- 
hatched region. Points below the horizontal hatched line are excluded by the requirement 
that M c + > 100 GeV. 
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Here Yf is the hypercharge of the field / and mj.f. = (f] + fj) 2 is the invariant mass of 
the fi, fj pair of fields. This matrix element can be simplified by assuming a common 
sfermion mass, Mj, and by putting in explicit values for the couplings 
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In the analysis of Ref. |I6|,P!7| the model for the heavy neutral lepton decay studied 



was based on a weak decay matrix element |24] 

M(N° - Vi £je£) = - 75)^.^7^(1 - l5 )u Vi . 



(10) 



If we compute the squared amplitude and average over the spin of the incoming heavy 
lepton we obtain 



\M\ 2 (N° -> Ui£j£+) = 16G 2 F (m 2 N0 + 
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So the distribution of the decay products from the R-parity violating decay will be exactly 
the same as the weak decay matrix element studied in []T5], |TT|| and therefore this model 



has exactly the same problem with the energy asymmetry Ae as that discussed in [ I6| , |l~7 



3.1 Neutralino Pair Production 

We simulated neutralino pair production using HERWIG 6.2 [p5| -|27|| .p| This allows us 
to simulate the production cross section with the correct momentum spectrum for the 
neutralinos and to determine whether they can decay within the NuTeV detector. Those 
events where the neutralino could decay in the detector were weighted with the probability 
that the neutralino decayed in the detector, for a given lifetime 

V^expl— , (12) 

where i = 1.4 km is the distance target-detector, Ax = 35 m is the length of the detector, 
(3c is the speed of the neutralino and r x o is its lifetime. The neutralino was further- 
more decayed with the full RPV matrix element 0,^]. We then applied the NuTeV 



kinematic cuts |T7j on the neutralino decay products. We required that the neutrali- 



nos decay within the fiducial volume^] (2.54 x 2.54 x 28) m 3 of the NuTeV detector at 

2 One modification to HERWIG was made in that we used the average of the central and higher gluon 
parton distribution functions from the leading-order fit of |||] . This will become the default in the next 
release of HERWIG. 

3 In the original version of our paper this number was smaller as found in (46). We thank T. Adams 
for drawing our attention to the corrected value in the published version [r^| . 
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Figure 5: Number of events in the NuTeV detector for neutralino production in B- meson 
decays as a function of the neutralino lifetime. 



a distance 1.4 km downstream of the production target. The muons produced in the 
neutralino decay were required to have energy, > 2.2 GeV, and the transverse mass, 
m T = \p T \ + ^/p2 + M 2 > 2 2 GeV, as in Ref. pHO]. Here P T and M T are the trans- 
verse momentum and mass of the visible decay products, respectively. 

As the production of a bino only occurs via t-channel squark exchange the cross section 
will depend on the (assumed degenerate) squark mass as ~ The number of events 

which would be observed in the NuTeV detector are given in Fig.^ as a function of 
the lifetime of the neutralino. Given the current limits on the squark mass from both 
30-031 and the Tevatron it is impossible, for any neutralino lifetime, to get 



LEP 

sufficient events to explain the NuTeV results via neutralino pair production. In Ref. |14 
the LEP constraints on the MSSM parameter space were not taken into account. 



3.2 Neutralino Production in 5-meson Decays 

As with the neutralino pair production we used HERWIG to simulate bb production. One 
of the produced B mesons was then forced at random to decay via RPV. The overall 
normalization was properly taken into account. The partial widths for the decays of the 
B and B + via RPV are given by 
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Figure 6: Regions in A232, A^ 13 parameter space in which we would expect 3 ± 1 events to 
be observed in the NuTeV detector. The limits from [0,0 on the couplings A' 113 (crosses) 



and A' 213 (diamonds) allow solutions between the two points for each of the masses shown. 
The region above the stars is ruled out for the coupling A' 213 by the limit on the product 
of the couplings A232A 213 from the limit on the branching ratio B + — > fi + u ||36|| . The 
hatched region shows the experimental bound on the coupling A^ 13 from perturbativity. 
The corresponding limits on the coupling A232 from both low energy experiments p,[T 
and perturbativity are not shown as our solutions do not extend into this region. 
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where p cm is the momentum of the decay products in the rest frame of the decaying meson, 
m u,d,b are the up, down and bottom quark masses respectively, m B o is the B° mass, ibb+ 
is the B + mass. Here Lf — —g'Yf L /2 for the left-handed fermions and Rf — g'Yf R /2 for 
the right-handed fermions. fs is the pseudo-scalar decay constant for B decays, M^o is 
the lightest neutralino mass, is the left down squark mass, M^ L is the left up squark 
mass, and is the right bottom squark mass. In Eqns. (0,0) we have assumed that 
the sfermions have a common mass Mj. The pseudo-scalar decay constant for the B 
system has not been measured experimentally and must be taken from lattice QCD. We 
have used the value 

f B = 204 ± 30 MeV, (15) 



from Ref. |J7| where we have added the errors in quadrature. The branching ratio for 
the decay B° — ► Xi^ was taken as an input and the branching ratio for B + — ► 
calculated from it using the above results. The same cuts were applied as in the previous 
section. The number of events which would be observed in the detector is shown in 
Fig.[|. This shows that even for branching ratios below 10~ 7 there is a significant range 
of neutralino lifetimes for which there are enough events to explain the NuTeV results. 
The present experimental upper limit on the branching ratio of the purely muonic decay 



is Br(B ± -f /z%) < 2.1 ■ 10" 5 p| 



Using the results for the RPV branching ratios of the B mesons and the neutralino 
lifetime we can find regions in (A232, A' 113 ) parameter space, for a given sfermion mass, in 
which there are 3 ± 1 events inside the NuTeV detector, this is shown in Fig. || We have 
included the low-energy bounds Eq.(^). In the case of the coupling A' 213 there is also a 
bound on the product of the couplings A232 • A' 213 from the limits on the branching ratios 
for B° — ► r~/z + and B + — > fi + u |[36|| , the latter giving the stricter bound 

A 2 2 13^3 2 />l+ ( 1 _ J^_\ 2 <2 Ax 1Q -B (16) 

32ttM| (m b + m u f V B+ \ m 2 B+ J 
Here Tb+ is the total width for the B + . This gives 

A^SS.SxlO-yy 2 . (17) 

In Fig. ^|we see that for every value of A^ 13 there are two solutions in A232, except for 
a minimum value of A^ 13 , below which there are no solutions. This can be understood as 
follows. The maximum fraction of neutralinos decays in the distant detector for a lifetime 
r = [3cr//£, i.e. when the decay length corresponds to the flight length, the distance 
between the production target and the detector. This optimised lifetime corresponds 
numerically to 

A, S9 = 5.3 • 10" 4 ( ~"i .. 1 . ( U 



( M f V 
V 200 GeV / 



This requires the minimum production rate and thus the minimum value of A- 13 , which 
is the dip in the curves in Fig. |6]. For larger values of A^ 13 the neutralino production is 
increased. We can then tune the lifetime of the neutralino such that the decay length 
is either shorter or longer than the flight length, yielding the two solutions shown in the 
figure. 



3.3 r-Decays 

As discussed in section |], in our model the neutralino can decay to jiru as well as /z/zza Us- 
ing the calculation of Eq.(|9]) we can compute the branching ratios 
Br m = Br(xi ^ I^I^^t) an d Br^ T = Br(xi ~^ T± H Tl> n)i which are displayed in Fig.[7|. 
For neutralino masses above 10-15 GeV the two decays have practically equal branching 
ratios. However, when the neutralino mass is close to the r-mass, Br^ T is phase space sup- 
pressed. For M x o = 5 GeV we have Br^ T = 0.287. In obtaining Fig. [7] the sfermions have 
been assumed to be degenerate and left/right stau mixing has been neglected.^ In princi- 

4 In models where the scalar masses are unified at the GUT scale the running of the masses to low 
scales forces the right stau to be lighter than the left stau. For low tan/? it is a good approximation to 
neglect left/right stau mixing. For large values of tan j3 the right stau becomes much lighter, but this does 
not contribute to the decay. It is thus a conservative assumption to require degenerate scalar fermion 
masses. 
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Figure 7: Branching ratios for the decay of a purely bino lightest neutralino via the RPV 
coupling A 2 32- The sfermions have been assumed to be degenerate and light /right stau 
mixing has been neglected. 



pie the NuTeV experiment can observe the (/ir z/)-modes through the decays: r — > evv 
and — > 7r ± (n • 7r°)z/, which would lead to unobserved (e,/i) and (vr,/i) events, respec- 
tively. Here (n ■ tt°) indicates an additional n = 0, 1, 2, 3 emitted neutral pions. Given the 
3 observed (fi, fi) events one would expect the following number of events for M x o = 5 GeV 



N (e,n) = 3-- §— ■ Br(r euu) w 0.21 (19) 

^) = 3 ■ ^ T ■ Er(r -> 7r(n ■ tt» « 0.56, (20) 

1 LJT urr 



where we have used the r branching ratios from Ref. [36]. Thus the non-observation 
of (e,/i)- and (71", \i) -events is consistent. We note that some of the r — ► 7r ± (n ■ 7r°)z/ 
decays would show extra activity in the detector and thus be rejected as pure 7r =1= events. 



Therefore the above estimate is conservative 38 



4 Future Tests of the R p Violating Model 
4.1 NOMAD 

The NOMAD experiment [j39| -f4l| was a neutrino oscillation experiment at CERN which 



was dismantled in 1999. The data however are still on tape and could be used to test the 
current proposal. We modified our program to estimate the event rate at NOMAD. For 
this we used the following numbers |3"9"|-pn| : distance target-detector: £ = 835 m, fiducial 
volume of the detector: V = (2.6 x 2.6 x 4) m 3 , target material: Beryllium, target density: 
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Figure 8: The predicted number of di-muon events at NOMAD as a function of the 
neutralino lifetime. We have used our model for neutralino production through 5-meson 
decays. We indicate the prediction for three different branching ratios of the neutral 
B-meson decay to neutralinos as in Fig.|5|. 

p = 1.85 g/cm 3 , target length: d = 1.1m, proton beam energy: E = 450 GeV, integrated 
number of protons: N p = 4.1 • 10 19 . Using these numbers we show our prediction for 
the number of events at NOMAD in Fig.[8]. For the same f?°-meson branching ratio we 
obtain about an order of magnitude more events than at NuTeV. Thus our model can be 
completely tested by the NOMAD data! 

The higher sensitivity at NOMAD can be understood as follows. The total bb produc- 
tion cross section for collisions on Be at NOMAD is 4.7 nb, while for collisions on BeO at 
NuTeV it is 94 nb. The total integrated luminosities are 5.58 • 10 11 nb -1 (NOMAD) and 
6.189 • 10 rib (NuTeV), respectively. Therefore the number of 66-events are 2.6 • 10 12 
(NOMAD) and 5.8 • 10 11 (NuTeV), respectively, i.e. about 4.5 times more at NOMAD. 
The NOMAD detector is closer and than the NuTeV detector and thus subtends a larger 
solid angle by about a factor of 3. The required neutralino lifetime is about the same be- 
cause NOMAD is about half the distance but the energy is also about half. The NOMAD 
detector is about 8 times shorter but the Lorentz boost is only about half the NuTeV 
boost, so this corresponds to a factor of 4. All in all we would expect about a factor of 
3.4 times more events at NOMAD than at NuTeV. Comparing Fig. [| with Fig. ||] we see 
that this is close to what the full numerical simulation gives. 

4.2 B-Factories 

As outlined above, for neutralino production we are relying on a rare B-meson decay 

B ± - ^ + Xi, (21) 



(22) 



which can possibly be observed at a present or future B-factory. In the Standard Model 
we have the decay — > / u ± + z/ with a predicted branching ratio |^| of about 3 • 10~'. 
This is probably just below visibility at BaBar |[42| . 

The decay (|2l|) differs from the Standard Model decay B ± — > xi + v in the energy of 



the //, which is now only £ M = (M|± + m\ - M*)/(2Mb±) « 0.27 GeV for M x = 5 GeV. 
We thus have a monochromatic muon with an order of magnitude less energy than in the 
Standard Model decay This is a distinctive signature which we propose for investigation 
at BaBar and other B-factories. We presume this is very difficult due to many sources of 
soft muons as background. Also the efficiency for such soft muons is typically very low, 
only about 5% 

The decay ( ^2j) is invisible, with the neutralino decay far outside the detector at a 
B-factory. If we have a B°-B° system and could tag one of the mesons, via a conventional 
decay, then we would have an unexpected invisible decay on the opposite side. We propose 
this as a possible signature for investigation by the experimental collaborations. 



5 Neutral Heavy Leptons 

In |16|, [FT] the NuTeV collaboration also considered the possibility of a neutral heavy 
lepton (NHL) to explain their observation. Here a NHL, Nn, i = 1, 2, 3 is considered as 
a primarily isosinglet field under SU(2)l with a small admixture of the light Standard 
Model neutrinos. This is discussed for example in Refs. f4"3"l,[44|l. We follow the notation 
of Ref. In general such a NHL has charged current (CC) and neutral current (NC) 



purely leptonic decays proceeding via a virtual W or Z -boson, respectively, 

N iL - £j+£t + u k , (CC) (23) 
N iL -> u m + £+ + £~. (NC) (24) 

For the NC-decay the charged leptons are from the same family, whereas for the CC- 
decay they can also be from different families. A given CC-leptonic decay is proportional 
to the mixing element [Uj^l 2 - There is a corresponding NC-decay proportional to the 
same mixing element for m = j. For a given set of NHL masses and mixings, we typically 
would expect both NC- and CC-decays to occur. k,n — 1,2,3 are free indices which all 
contribute to the decay rate, independent of the mixings. 

NuTeV observe an excess of di-muon events. Assume we have one NHL, iVj, with mass 
M^. = 5 GeV, and the other NHL's unobservably heavy. The di-muon events could occur 
through CC-decays with j = k = 2 and the mixing element L^TVj or through NC-decays 
with n = 2 and the mixing elements UmNi, m = 1,2,3. For j = k = n = 2 we obtain 
di-muon events through both NC- and CC-decays. 

If the CC-decays contribute, i.e. j = 2, we would expect there to be accompanying 
(e, fi) events with similar probability, from the k = 1 mode. For example for a non- 
vanishing \U2Ni \ 2 , using the decay rates given in pHfl , we obtain the ratio of (e, fi) to (//, /x) 
events given by R e ^/^ = T(A r 2 — ► e + /i _ ^ M )/r(A r 2 — > /x 4 ' We plot this as a function 
of the NHL-mass in Fig.|9|. From the plot we see that we would expect more (e/i) events 
than (/i, /x)-events. This is excluded by the NuTeV non-observation of such events. 

If the NC-decays contribute we can expect further (e, e) and (r, r) events. The latter 
are kinematically suppressed as in the fi p case above. A search for the (e, e)-modes has 
only been presented for low-energy electrons . However, given a non- vanishing mixing 
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Figure 9: Ratio of (e, //)- to (//, /i)-e vents, -R, 



A^, versus the mass of iVj 



e/Vw — r(Af 2 ->Ai+/i-^ M ) 



in the decay of a NHL, 



element U m Ni which gives the (//, /i)-events via NC-decays we would expect further CC- 
decays: Nn — > if + i\ + Vk-, k = 1,2,3. In particular, for j = 1,2 this leads again to 
(e, /x) events which were not observed. 

We have thus eliminated all cases except a special model, which we consider in more 
detail. Assume j — 3 and is the only non-negligible mixing element. Furthermore, 
as above, assume Mat. = 5 GeV and the other NHL's are unobservably heavy. We then 
have the following decay modes 

N {TTUr^fiu^reUe} , (25) 
Ni {u T ee, z/,-/J/i, v t tt} ■ (26) 

The r and rr decay modes are kinematically suppressed as in the ft p case discussed above 
and the observed di-muon events are obtained form the NC-decay. This model has been 
studied by the NOMAD collaboration for M N . = 10 - 190 MeV |TJ. 

We now estimate the event rate for this model (j = 3, UsN t 7^ 0). The production 
mechanism will go either via (CC or NC) Drell-Yan with the tau neutrino mixing with the 
Ni or via B-meson decays. We have computed the Drell-Yan production cross section to 
be gdy = (9(10 _1 pb) • | f^3iVi 1 2 - The neutral current contribution to the NHL production 
is more than an order of magnitude smaller. The total integrated luminosity at NuTeV 
corresponds to ~ 6.2 • 10 6 pb _1 giving the number of produced Nf N^° d ~ 4 • 10 6 \U 3 Ni\ 2 - 
Of these only about 1% fly in the direction of the detector JTJJ, leaving us with ~ 



4 • 10 4 It/ajvJ 2 . In order to estimate the total number of events we must combine this with 
the fraction of Aj which decay in the detector given by Eq. (|T2"D. The total event rate is 
proportional to 

N ev « N Ni exp {- a \U 3N f} ■ b\U 3N f = 4 • 10 4 6 \U 3Ni | 4 exp {- a \U 3Ni \ 2 } (27) 

where a = £ / ' (f3'~ic{TN i \U3N i \ 2 ) and b = Ax / (PjcTNi \U 3 Ni | 2 ) from Eq. (|12^) are independent 
of I^atJ. The event rate is maximal for 1 1/3 ^ | 2 = 2/ a. We obtain an upper limit on the 



lifetime if we assume the NC decay is dominant. The latter we determine through the 
scaled muon lifetime 



r Ni < r, \U 3Ni \- 2 = 9 • 10~ 15 s \U 3Ni \~ 2 . (28) 

We then obtain a = 5.2 • 10 8 /7 and b = 1.3 • 10 7 /7. For 7 = 10, for example, we obtain the 
maximal event rate for {U^^ | = 9 • 10 -5 , which is compatible with the independent bound 
|t^3JVi| 2 < 0.016 (45|. Following Eq. fll2]) the total fraction decaying in the detector is 
then roughly 1.1%. Combining this with the previous estimate of the number produced 
we get a total maximal number of events of about N™ ax = 5 ■ 10~ 7 , which is of course too 
small. 

The reason this is so much smaller than in the supersymmetric model is that the 
lifetime of the NHL is typically much shorter. Thus the NHL's typically would decay 
well before the detector. We get the maximal number of events when the lifetime is 
approximately the flight time. For this we need a very small I^atJ. Since we only have 
one parameter in this model this feeds into the cross section resulting in the very low rate. 
We do not expect the production via 5-mesons to help. The branching ratio is suppressed 
compared to the SM decay branching ratio Br(B + — > r + Ni) = \U^Ni\ 2 Br(B + — > t + v t ) w 
7 ■ 10 _5 |£/3Ar I 2 and thus also too small. 



6 Conclusions 

We have reconsidered the NuTeV di-muon observation in the light of supersymmetry 
with broken R-parity and neutral heavy leptons. We have shown that it is not possible 
to obtain the observed event rate with pair production of light neutralinos or via the 
production of neutral heavy leptons. However, we have introduced a new production 
method of neutralinos via 5-mesons. Due to the copious production of i?-mesons in the 
fixed target collisions the observed di-muon event rate can be easily obtained for allowed 
values of the R-parity violating couplings. 

The model we have proposed can be completely tested using current NOMAD data. 
We suspect this is true of any model one might propose. If the NOMAD search is negative 
our model is ruled out and the NuTeV observation is most likely not due to physics beyond 
the Standard Model. 

It is worth pointing out that through this mechanism we have opened a new sensitivity 
range in the R-parity violating couplings. At colliders we can probe the range where the 
neutralino decays in the detector. For a photino neutralino this corresponds to 

r ^( rh V/i00GeV\ 5/2 / fa \ 2 f5GeV\ 5/2 

(29) 

Here we have substituted the light neutralino mass we are considering. For significant 
boost factors we thus can probe couplings at most down to 1CT 3 . From Fig. || we see that 
for a 200 GeV sfermion we can probe couplings down to about 5 • 10~ 6 , which is more 
than two orders of magnitude smaller! It is thus worthwhile to study the production of 
neutralinos via mesons at fixed target experiments in more detail. 

Before concluding we also note that one might worry that the lightest supersymmetric 
Higgs boson would decay dominantly to the two light neutralinos and thus be invisible. 
However, as with the Z° boson, the Higgs does not couple to a Bino neutralino. 
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